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ABSTRACT 

Luminous Red Galaxies (LRG) from the Sloan Digital Sky Survey are considered among the 
best understood samples of galaxies, and they are employed in a broad range of cosmologi¬ 
cal studies. Because they form a relatively homogeneous population, with high stellar masses 
and red colors, they are expected to occupy halos in a relatively simple way. In this paper, 
we study how LRGs occupy massive halos via direct counts in clusters and we reveal several 
unexpected trends suggesting that the connection between LRGs and dark matter halos may 
not be straightforward. Using the redMaPPer cluster catalog, we derive the central occupa¬ 
tion of LRGs as a function richness, iVcen(A). Assuming no correlation between cluster mass 
and central galaxy luminosity at fixed richness, we show that clusters contain a significantly 
lower fraction of central LRGs than predicted from the two-point correlation function. At halo 
masses of we find A^^en = 0.73, compared to A^^en of 0.89 from correlation studies. 

Our central occupation function for LRGs converges to 0.95 at large halo masses. A strong 
anti-correlation between central luminosity and cluster mass at fixed richness is required to 
reconcile our results with those based on clustering studies. We also derive Pbnc^ the proba¬ 
bility that the brightest cluster member is not the central galaxy. We find Pbnc ~ 20 — 30% 
which is a factor of ~ 2 lower than the value found by |Skibba et al. ( 201 l| l. Finally, we study 
the radial offsets of bright non-central LRGs from cluster centers and show that bright non¬ 
central LRGs follow a different radial distribution compared to red cluster members, which 
follow a Navarro-Frank-White profile. This work demonstrates that even the most massive 
clusters do not always have an LRG at the center, and that the brightest galaxy in a cluster is 
not always the central galaxy. 


Key words: galaxies: clusters: general 


1 INTRODUCTION 


Luminous Red Galaxies (LRGs) are early-type massive galaxies 
consisting mainly of old stars with little ongoing star formation. 
Because LRGs are very luminous and can be detected to z ~ 0.5, 
they are commonly used for studies of large scale structure. The 
Sloan Digital Sky Survey (SPSS, [York et al.|2000^ targeted a large 
sample of LRGs out to z ~ 0.5 using color and magnitude cuts 
described in |Eisenstein et al.1j200l^ . LRGs from SDSS have been 
used for a wide variety of purposes, such as the detection of baryon 
acoustic oscillations (e.g., Eisenstein et al.|2005 Kazin et al. 2010| (, 
clustering studies (e.g.,|Zehavi et al.|200h Reid & Sperge |2009[ 


jPadmanabhan et al.|2009[ Parejko et al.|2013|, weak tensing and 


cross-correlation studies (e.g., [Mandelbaum et al.||2006 

IHikage 

et al.|2013 

I, and redshift space distortion (RSD) studies i 

Cabre & 

Gaztanaga 

2009 [Reid et al.|2009||Samushia et al.|2012[|Hikage & 

Yamamoto 

2013b. 


All of the studies described above require an understanding of 
how LRGs trace dark matter halos, and, therefore, large scale struc¬ 
ture. One important aspect of how LRGs trace dark matter is how 
often they coincide with the centers of dark matter halos. The iden¬ 
tification of central galaxies in cluster samples is critical for a va¬ 
riety of studies, including the determination of cluster halo masses 
with weak gravitational leasing ( [George et al.|2012l ), and RSD stud¬ 
ies (jHikage & Yamamoto|2013)). Because LRGs are luminous and 
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trace massive overdensities, it is often assumed that central galax¬ 
ies are also LRGs. Indeed, in most galaxy clusters with at least 
one LRG, the Brightest Cluster Galaxy (BCG) is an LRG. Further¬ 
more, it is usually assumed that BCGs are central galaxies ([van den| 
|Bosch et al.|2004[|Weinmann et al.|2006[[Budzynski et al.|2012| l. 

However, recent results suggest that the central galaxy is not the 
BCG in all clusters ( |van den Bosch et al.|2005[|Coziol et al. 2009 


Sanderson et al.|2009||Skibba et al.|201 l||Einasto et al.|201 1| 2012 


Hikage et al.|2013l Sehgal et al.|2013[|Lauer et al.|2014^ . For exam¬ 
ple, by analyzing differences between the velocities and positions 
of BCGs relative to other cluster members, [Skibba et al.| j201 1^ 
found that 40% of BCGs may in fact be satellite galaxies. 

Knowing how LRGs occupy dark matter halos, and especially 
whether or not the most dominant LRGs coincides with cluster cen¬ 
ters, is of particular importance for RSD studies. |Hikage (& YaJ 
|mamoto| j2013| l found that satellite LRGs have a significant effect 
on the higher order multipoles of the LRG power spectrum. This 
is simply due to the fact that satellite galaxies have large peculiar 
velocities that cause structures to be smeared out along the line-of- 
sight (the “Finger-of-God” effect). Hence, understanding if LRGs 
trace halo centers and how often LRGs are satellites is an important 
ingredient for large scale structure studies with LRG samples. 

Halo Occupation Distribution models (HOD: Seljak J2000) 


Peacock & SmithPOOPlISi 


et al.|20()Tl|Berlind & Wein-| 


coccimarro 


berg 2002^ are the most popular framework for describing how 
LRGs populate dark matter halos. The HOD model relies on the 
choice of a functional form for the probability that a halo of mass 
Mh contains N objects of a particular type (P{N\Mii)). This 
probability is typically divided into a contribution from central 
and from satellite galaxies. The mean probability for each com¬ 
ponent is the central occupation function, A'cen, and the satellite 
occupation function, Nsm, respectively. At group and cluster scales 
(Mh > lO^^^M©), the central occupation function of LRGs is often 
assumed to converge to unity. If this assumption is correct, then the 
central galaxies of massive halos should all be LRGs. One of the 
primary goals of this paper is to test this assumption. 

HOD models for LRGs are typically constrained by mea¬ 
surements of the LRG two-point correlation function (e.g., |Zheng| 
|et al.|20()5]|Masjedi et al.|2006^ or by u sing the counts-in-cylinders 
method developed by Reid & Spergel|p009| l. However, a comple¬ 
mentary way to measure the HOD of LRGs is to directly count 
LRGs in clusters as a function of halo mass. For example, |Ho| 
|et al.| j2009^ studied the HOD of LRGs using a sample of 47 
clusters from the ROentgen SATellite (ROSAT) X-ray survey. The 
main caveat to this direct approach is that it relies on the avail¬ 
ability of a trustworthy cluster catalog with known centers and 
halo masses. Early cluster catalogs from SDSS such as MaxBCG 
(|Koester^e£^aLp007j(^are known to have problems with centering. 
[Johnston et al.| j2007| > estimate that the mis-identification of cen¬ 
trals in the MaxBCG catalog is roughly ~ 30%. Recently, how¬ 
ever, much progress has been made in understanding how to im¬ 
prove centering algorithms for cluster finders ([Rykoff et al.|2014| 
[Rozo & Rykoff|2014^ . 


The goal of this paper is to compare the central occupation 
function of LRGs as inferred from galaxy correlation measure¬ 
ments to direct counts of central LRGs in clusters. We use the 
redMaPPer cluster catalog ( [Rykoff et al.||2014l [Rozo & Rykoff[ 
[2014[ l, which includes carefully selected central galaxies based on 
luminosity, color and local galaxy density. One of the advantages 
of redMaPPer over previous cluster finders is that it adopts a proba¬ 
bilistic approach to cluster centering. This probabilistic approach is 
especially useful when, for example, two clusters are merging and 


when the central galaxy is not obvious. We caution that while the 
comparison of X-ray to optical centers in galaxy clusters so far sug¬ 
gests that the redMaPPer centering probabilities are accurate ( [Rozo[ 
[& Rykoff|[2014^ , the samples used in these comparisons are still 
small, and rely heavily on X-ray selected sub-samples of clusters. 
More detailed investigation into the redMaPPer centering probabil¬ 
ities is still warranted, but for the purposes of this study we will 
assume the centering probabilities from the redMaPPer algorithm 
are correct. In Paper 11, we perform additional tests on the redMaP¬ 
Per centering probabilities by using weak gravitational lensing and 
cross-correlations. Preliminary results from this work indicate that 
the redMaPPer centroids are, on average, better indicators of halo 
centers than BCGs, adding confidence to the redMaPPer centering 
probabilities. 

This paper is organized as follows. In Sj^we briefly describe 
the redMaPPer cluster catalog and our photometric SDSS LRG se¬ 
lection. In we present the central occupation function of LRGs 
as inferred from the redMaPPer cluster catalog, and compare four 
different definitions of the central galaxy. Our main results are pre¬ 
sented in O Finally, we discuss the results and draw our conclu¬ 
sions in ^4]and 

Throughout this paper, we adopt a flat ACDM cosmology with 
h = 0.7, flm ~ 0.25, and erg = 0.8. This cosmology is chosen 
for consistency with the redMaPPer richness-halo mass calibration 
in [Rykoff et al.[ j2012^ . All distances are expressed in physical 
Mpe units. Halo mass is defined as M 2 oob = M(< J? 200 b) = 
200p|7ri?2oob> where J? 200 b is the radius at which the mean inte¬ 
rior density is equal to 200 times the mean matter density (p). All 
magnitudes are given on the AB system. 


2 DATA 

2.1 RedMaPPer Cluster Catalog 


The red-sequence Matched-filter Probabilistic Percolation 
(redMaPPer) algorithm is a red-sequence based cluster-finder 
optimized for identifying clusters in large multi-wavelength 
surveys. The details of the algorithm can be found in [Rykoff] 
[et al.[ ( [2014[ ( and [Rozo & Rykoff[ ( [2014[ (. Comparisons between 
redMaPPer and other cluster finding algorithms are performed 
in [Rozo & Rykoff[ ( [2014) ; [Rozo et aL]p014[ l. The key elements 
of the redMaPPer algorithm are briefly summarized here. Using 
galaxies with spectroscopic redshifts, redMaPPer first constructs 
a red-sequence model as a function of luminosity and redshift. 
Once the red-sequence model is calibrated, redMaPPer then uses 
an iterative method to identify clusters assuming simple radial and 
luminosity filters. The redMaPPer algorithm is fully probabilistic. 
All potential cluster members are assigned a probability, Pmem, of 
being a red cluster member. Cluster membership probabilities are 
tested using spectroscopic redshifts and are accurate to 1% ( [Rozo] 
[et al.[|20T4) . Each cluster is assigned a photometric redshift, de¬ 
noted zx, that is estimated from high-probability cluster members. 
Finally, redMaPPer estimates cluster richness (denoted A) as the 
sum of membership probabilities over galaxies: 




Pfree,i Pmem,i 


i 


( 1 ) 


where pfee represents the probability that a galaxy is not a member 
of another cluster (in general, p^ee ~ 1). The sum is performed 
over all galaxies within a characteristic cutoff radius, Rx, which is 
a richness-dependent aperture that minimizes scatter in the mass¬ 
richness relation. 
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Luminous Red Galaxies in Clusters 3 


A key feature of the redMaPPer algorithm is that cluster cen¬ 
tering is also probahilistic. For each cluster, redMaPPer provides a 
list of the five most likely central galaxies. Each of these five cen¬ 
tral galaxy candidates is assigned a probability, Pcen, of being the 
cluster center. Centering probabilities are defined using three fil¬ 
ters: a luminosity filter (i^cen), a photometric redshift filter (Gcen), 
and a local galaxy density filter (/ccn). The product of these three 
terms determines the overall centering filter: 

t^cen ^cen (rrlilzx, A)Gcen(2red)/cen(w|2A, A), (2) 


closer than 55 arcsec (SDSS-I/II) or 62 arcsec (SDSS-III). Large 
clusters may contain several LRGs and are hence more likely to be 
affected by fiber collisions than smaller systems, which have only 
one LRG. As a result, in the redMaPPer catalog, 23% of LRG clus¬ 
ter members lack a spectroscopic redshift. In order to circumvent 
this issue, we will use photometrically-defined LRG samples. We 
use photometric data from SDSS DR8 for both the classical LRG 
and LOWZ selections. 

In this paper, we consider the following three LRG samples: 


Here, rrii is the i-band apparent magnitude, Zid is the redshift of 
an individual galaxy (inferred from the color), w is the local galaxy 
density. An important point to note is that the redshift filter Gccn 
is broader than the cluster red-sequence filter and allows galaxies 
with slightly offset colors, and, therefore, lower Pmcm values, to be 
considered central galaxy candidates. This allows for the possibility 
of residual amounts of star-formation, and, therefore, bluer colors, 
in centrals. For each cluster, centering probabilities are normalized 
as: 

5 

X] Pcen.i = 1 (3) 

i=l 

One consequence of Equation is that the brightest cluster 
member is not necessarily the galaxy with the highest value of Pcen- 
This is because, in addition to luminosity, the central galaxy selec¬ 
tion also considers local galaxy density via the fcm filter. Hence, 
in some cases, a less luminous galaxy may have a higher centering 
probability because it is spatially coincident with the dense cluster 
core. 

We use the redMaPPer v5.10 cluster c atalog jR ozo et a l.|2014^ 
based on the SDSS 8th Data Release (DR8: [Mhara et al.|2011t . We 
select clusters in the range A > 20 and 0.16 < 2 < 0.33. The 
lower redshift limit, 2min = 0.16, is set to avoid low redshifts where 
the LRG color cuts begin to select a fainter population of galaxies 
( [Eisenstein et al.|2001| l. The upper redshift limit, 2max = 0.33, is 
set so that the resulting sample of clusters is volume limited. In 
addition, we select clusters for which fewer than 20% of the cluster 
is masked (raask_f rac <0.2). Altogether, there are 7,730 clusters 
in this redshift and richness range. 

Our SDSS DR8 redMaPPer cluster catalog is limited to i < 
21.0. This is roughly one magnitude deeper than LRG samples 
(see Section [23^ , ensuring that our cluster membership selection 
is complete for LRG-type galaxies. 


(i) DR8PhotLRG : photometrically defined “classical” LRGs 
selected from DR8 photometry. 

(ii) DR8PhotLRG-MgCut: same as above but with an additional 
g-band absolute magnitude cut (see Section [Z3) . We use this sam¬ 
ple to compare with previous results derived from clustering stud¬ 
ies. 

(iii) PhotLOWZ : photometrically defined BOSS LOWZ sample 
selected from DR8 photometry. 


We base our LRG selections on DR8 photometry but we note 
that sub-sets of the LRG samples were targeted based on photom¬ 
etry from earlier data releases. In Appendix A, we compare our 
samples based on DR8 photometry to LRG samples based on DR7 
photometry and to the original spectroscopic target catalogs. In ad¬ 
dition to covering a larger area than DR7, the DR8 photometry uses 
a new sky subtraction algorithm, affecting mainly bright galaxies 
(r < 16) ([Abazaj ian et al.|2009^ , and a new calibration ( [Padman^ 
|abhan et^ 2008| l. While the total numbers of LRGs selected using 
DR7 and DR8 photometry are the same, the actual galaxy samples 
differ at the ~ 5% level. This is due to scatter around the color 
and magnitude cuts imposed by the LRG selections below. As we 
show in the Appendix A, the differences in photometry for LRGs 
are small (clipped RMS Snir = 0.07) and there are no systematic 
offsets in magnitudes or color. The lack of systematic offsets sug¬ 
gests that the improved sky subtraction in DR8 has little effect on 
LRGs. This is expected since most LRGs are fainter than r < 16, 
where the effects of improved sky subtraction are most noticeable. 
Therefore, while the DR7-selected and DR8-selected LRG samples 
are not identical, they are comparable. 


2.2 Luminous Red Galaxies 

LRGs are intrinsically red and bright galaxies selected from SDSS 
( [Eisenstein et al.|[2001 ^. LRG samples extend fainter and farther 
than the SDSS main galaxy sample ( [Strauss et al.|200^ and con¬ 
sist of mainly passively-evolving, massive, early-type galaxies with 
redshifts in the range 0.16 < 2 < 0.7 (though we restrict the red¬ 
shift range to 0.16 < z < 0.33 here). In this paper, we will distin¬ 
guish “classical LRG” samples, which were targeted as part of the 
SDSS-I and SDSS-II programs ( [Eisenstein et al.[2001[ (, from LRGs 
that were targeted as part of the SDSS-III Baryon Oscillation Spec¬ 
troscopic Survey (BOSS) program ( [Dawson et al.|20T3) l. We will 
refer to LRGs selected in SDSS-I and II as “Classical LRGs”. We 
will use the term “LOWZ” to refer to low redshift (z < 0.4) LRGs 
from the SDSS-III BOSS program. 

Fiber collisions occur because two fibers on the same SDSS 
spectroscopic plug plate (referred to as a “tile”) cannot be placed 


2.3 Classical LRGs and DRSPhotLRG Cluster Member 
Selection 

In this section we explain how we construct our sample of classical 
LRG cluster members, the “DRSPhotLRG” sample. As stated in 
Section 113 we only consider cluster members with cluster red¬ 
shifts between 0.16 < 2 < 0.33. This redshift range is simi¬ 
lar to the range considered for a variety of LRG studies: [Zehavi[ 
[et al.[ ( [2005[ l use 0.16 < 2 < 0.44, Reid & Spergel[ f2009 1 use 
0.16 < 2 < 0.36, and [Kazin et al.[j2010| use 0.16 < z < 0.36. 

Classical LRGs are selected using a series of color-color and 
luminosity cuts in SDSS. Because the 4000A break moves from 
g-band to r-band at 2 ~ 0.4, for classical LRGs there are two sep¬ 
arate selections at redshifts above and below 2 ~ 0.4. The LRG 
selection at 2 < 0.4 is called “Cut I”. Given our redshift range, 
in this paper, we only consider “Cut I” LRGs. All colors are de- 
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fined using model magnitudeQand all quantities are corrected for 
galactic extinction jSchlegel et al.|19^ . 

We query the SDSS PhotoPrimary table (see Appendix B 
for the full query) for all cluster members to obtain the photomet¬ 
ric quantities necessary to apply the LRG “Cut I” selection. Here¬ 
after, the subscript “mod” denotes model magnitudes, which are 
derived by adopting the better fitting luminosity profile between 
a de Vaucouleurs and an exponential luminosity profile in the r- 
band. The subscript “cmod” denotes composite model magnitudes, 
which are calculated from the best-fitting linear combination of a 
de Vaucouleurs and an exponential luminosity profile. For classi¬ 
cal LRGs, luminosity cuts are defined using rpetro, the Petrosian 
r-band magnitude. Two ancillary colors, c±_ and cy, are defined to 
align perpendicular and parallel to the locus of LRG galaxies in the 
(^Kiod rmod, Pmod fmod) plane! 


C_L — (Pmod imod) (f?mod r’mod)/4 0.177, (4) 

C|| = 0.7(gmod - fmod) + 1.2[(rmod “ imod) “ 0.177]. (5) 


Galaxies at 2 < 0.4 are distributed close to the linear locus c± = 0, 
and C|| denotes where galaxies fall along this locus. Galaxies at 
higher redshifts have larger cy. 

In this paper, we focus only on LRG samples at 2 < 0.4 which 
are selected via “Cut I”, as defined on the DR8 target selection 
webpag^ 


^’Petro 

< 

13.116-fcy/ 0 . 3 , 

(6) 

^’Petro 

< 

19.2, 

(7) 

|cx| 

< 

0.2, 

(8) 

/^50 

< 

24.2 mag arcsec”^, 

(9) 

'^mod 


0 

to 

(10) 


where /rso is the average surface brightness inside a Petrosian half 
light radius, and rpsf is the point spread function (PSF) magnitude. 
Equation l[§ ensures that the LRG sample has a roughly constant 
absolute magnitude limit between 0.16 < 2 < 0.4. Equation 0 is 
set by the S/N requirement for spectroscopic redshift measure¬ 
ments using the SDSS telescope. The c± cut limits the redshift 
range to 2 < 0.4. The cut removes very low surface bright¬ 
ness objects, which may have odd colors. The final cut is a star- 
galaxy separator. As shown in Table ^ Equations 0 and 0 have 
the most impact in defining the LRG cluster member sample. How¬ 
ever, equation 0 only has a small impact on the LRG central sam¬ 
ple simply because centrals are typically brighter than this limit. 

With the extinction-corrected composite magnitudes, we cal¬ 
culate absolute magnitudes Mugriz for cluster members as follows: 

Mugriz = mugriz — [k-coiTection] — [distance modulus]. 


We use the software package k-correct version 3.1 to perform 
this calculation jBlanton & Roweis|2007[ >. To calculate the distance 
modulus, we use a spectroscopic redshift when it is available, oth¬ 
erwise we use the cluster redshift, 2 a . Absolute magnitudes derived 
with this equation are used in Eigures[^[^and[5 

Figure [^illustrates how Equations Q andjVj operate. In this 


^ Note that model magnitude were improved in SDSS DR2 and later data 
releases. This caused changes in galaxy colors (gmod — fmod) and (r^od — 
imod) of about 0.005 magnitude. The DRl LRG selection criteria originally 
specified by |Eisenstein et ar]j2001j would lead to an increase in the LRG 
number density by about 10% using DR2 photometry. To account for this, 
the LRG selection was slightly modified. 

^ http://www.sdss3.org/dr8/algorithms/target_selection.php 


Table 1. Percentage of cluster members and central galaxies that fail to 
pass the LRG selection cuts given by Equations §-(Tg. For example, 
among all 507,874 cluster member galaxies, 492,949 (97.1%) do not sat¬ 
isfy Equation and 390,840 (77.0%) do not satisfy Equation j^. Each 
cut is applied separately in this Table. Applying all of these cuts 12,001 
cluster members (2.36%) and 5,823 central galaxies (75.3%) are classified 
as LRGs. 


LRG cut 

Eq.0 

Eq.0 

Eq.jS) 

Eq.0 

Eq.® 

All cluster 
members 

97.1 

77.0 

5.37 

1.65 

2.04 

Most likely 
central galaxies 

24.4 

0.776 

0.285 

0.181 

0.0129 


figure, we show how cluster members and central galaxies from 
our sample populate the rpetro versus cy plane. There are additional 
cuts that define LRGs (Equations [TO^ but in practice these ad¬ 
ditional cuts only play a minor role (see Table [T](. The left panel 
in Eigure[T]shows cluster members color-coded by absolute magni¬ 
tude and illustrates that Equation 0 (diagonal line) roughly selects 
a population of LRGs brighter than an absolute magnitude limit of 
Mi ~ —23.2. The right hand side of Eigure^ shows redMaPPer 
centrals (the most likely central galaxies according to redMaPPer). 
Roughly 25% of redMaPPer centrals do not satisfy Equations 0- 
| |1Q[ (, mainly because they are fainter than Mi ~ —23.2. 

In order to obtain a volume-limited LRG sample, clustering 
studies often apply an additional cut based on the gr-band absolute 
magnitude . Mg jReid & Spergel|2009[ |Zehavi et al.|2005[|Kazin| 


etal. 


2010 1 . After this Mg cut, the LRG number densit^i 


Its riLRG = 


3.4 X 10 ® Mpc To make our sample consistent with 


Reid & 


|Spergel| ( |2009[ l, we apply a Mg cut of —23.2 < Mg < —21.2 
to the “DRSPhotLRG” sample to construct the “DRSPhotLRG- 
Mgcut” sample. E or consistency, we fo llow the same method to 
compute Mjjas in ] 


Kazin et al. 


l 2010|- 


Mg = rpetro—[distance modulus] — Ag 

4“ (0 ) 2:calibration ■ 


( 11 ) 


Since the observed r-band is close to the rest-frame ^-band, 
rpetro can be used to predict Mg. Ag and {g — r) are color and k- 
corrections from Table 1 in |Eisenstein et al.| ( |2001^ , and 2caiibration = 
0.2 accounts for evolution in the rest-frame colors from 2 ~ 0.3 to 


rived from Kazin et ^ ( |2010[ l yield similar results, despite the fact 
that |Kazin et al. ( 20 10^ used DR7 photometry and we use DR8. We 
find the difference in Mg is small (RMS SMg = 0.07), consistent 
with the scatter in the photometry and that the differences in Mg 
do not affect the number of LRGs. Therefore, our results can be 
compared to results that directly use the [Kazin et al.| ( |2010^ LRG 
selection. 


® Our values assume h = 0.7. In h inverse units, the number density of 
cla.ssical LRGs is npRo = 1-0 X 10“^ (/i“^Mpc)“® 

When quoting values for M; we assume h = 0.7. However, our Mg 
values assume h inverse units for consistency withjKazin et al.|^2010J and 
[Reid & Spergel|(|2009f 

^ http://cosmo.nyu.edu/~eak306/SDSS-LRG.html 
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2.51 . 

■ Cluster members 


2.0 - 


+ M,<-24.14 



+ -24.14 < Mi <-23.79 
+ -23.79 < Mi <-23.43 
+ -23.43 < Mi <-23.07 


o 

Di 

-23.07 <Mi<-22.71 
+ -22.71 <Mi<-22.36 
+ -22.36 < Mi <-22.0 




redMaPPer central galaxies 


+ LRG centrals 
+ non LRG centrals 


o 



^Petro 


^Petro 


Figure 1. Left panel: C|| as a function of Petrosian r-band magnitude, rpetro? for cluster members with 0.16 < 2 : < 0.33, Mi < —22 and A > 20 (data points 
are randomly down-sampled for visual clarity). Cluster members are color coded by their 2 -band absolute magnitude, M^. The classical LRG cuts select cluster 
members brighter than Mi ~ —23.2. Among all cluster members, about 2% are classified as LRGs. Right panel: cy versus rpetro for the most likely central 
galaxies chosen by redMaPPer. Red points show redMaPPer central galaxies that are LRGs while black points show centrals that fail the LRG selection. Most 
of the central galaxies in our sample that fail to pass the LRG cuts are intrinsically fainter than Mi ~ —23.2 (also see Figure]^. Roughly 25% of central 
galaxies in our cluster sample fail the LRG cuts. Central galaxies are primarily rejected from the LRG sample by Equation j^. 


2.4 BOSS LOWZ Photometric Cluster Member Selection 

In this section we explain how we construct our sample of BOSS 
“LOWZ” photometric cluster members (the “PhotLOWZ” sample). 
The details of the LOWZ selection can be found in IDawson et al.l 
( [2013^ . The main difference between classical LRGs and LOWZs 
is that the BOSS selection extends about 0.4 magnitudes deeper 
than the classical LRG selection. As a result, the LOWZ sample 
has a higher number densit)[^than classical LRGs : ?^LOwz = 1-0 X 
10“^ Mpc“®. In addition, the LOWZ selection includes a bright 
magnitude cut at rcmod > 16 which did not exist in the classical 
LRG selection. Finally, the LOWZ selection uses composite model 
magnitudes for magnitude cuts instead of Petrosian magnitudes. 

Again, to avoid issues with fiber collisions, we define a pho¬ 
tometric LOWZ sample. The ancillary colors c±,lowz and cy^Lowz 
are different than the classical LRG selection criteria: 

C_L,LOWZ — (fmod tmod) (jlmod r’niod)/4 0.18, (12) 

t-||,LOWZ ~ 0.7((^niod t'mod) A 1.2[(7‘niod fmod) 0.18]. (13) 

With these new definitions we select PhotLOWZ galaxies using 
these cuts: 

rcmod < 13.5 + C||/0.3, (14) 

16 < rcmod < 19.6, (15) 

|ci|<0.2, (16) 

rpsf - rcmod > 0.3. (17) 

As in the classical LRG selection, the c± cut crudely selects galax¬ 

ies in the redshift range 2 < 0.4, Equations and © give an 
absolute magnitude and an apparent magnitude cut, respectively, 

® Our values assume h = 0.7. In h inverse units, the LOWZ number den¬ 
sity is T2LOWZ = 3.0 X 10“^ (h~^Mpc)~^ 


Table 2. Percentage of cluster members and central galaxies that fail to pass 
the LOWZ selection cuts given by Equations (H-Oz}- In total, before any 
cuts, there are 507874 member galaxies and 7730 central galaxies. Among 
all the cluster members, 479,040 (94.3%) do not satisfy Equation 0 and 
316,028 (62.2%) do not satisfy Equation 0. Applying all of these cuts, 
24,761 cluster members (4.88%) and 7,001 central galaxies (90.6%) are 
classified as LOWZs. 


LOWZ cut 

Eq.0 

Eq.0 

Eq. (ig 

Eq.0 

All cluster 
members 

94.3 

62.2 

5.37 

4.74 

Most likely 
central galaxies 

7.84 

1.54 

0.285 

0.0259 


and Equation l |17[ > is the star-galaxy separator. As shown in Table 
1^ Equation ( |14^ is the most stringent cut in defining the PhotLOWZ 
cluster member sample. 

In Eigurel^we illustrate how the LOWZ selection operates. 
As in Eigure [^cluster members color-coded by absolute z-band 
magnitude are shown on the left and central galaxies are shown 
on the right. Absolute magnitudes are calculated in the same way 
as in Section [23l As can be seen by comparing Eigures[2and|^ 
the LOWZ selection reaches ~ 0.3 magnitudes fainter in the i- 
band than the classical LRG selection. At the bright end, 1.3% (one 
hundred galaxies) of the most likely redMaPPer central galaxies are 
rejected by the rcmod > 16 cut. 

2.5 Selecting LRGs in clusters 

Eor most of our calculations, we consider all LRGs in a clus¬ 
ter with nonzero membership probability, which we will call the 
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Figure 2. Same as Figure[^but for the BOSS LOWZ sample. Compared to classical LRGs, LOWZ galaxies extends to a fainter limit of Mi ~ —22.8. The 
right panel shows that the LOWZ selection includes a bright cut at rcmod > 16, which excludes 1.3% of the central galaxies in our sample. 


■ Cluster members 

+ M,<-24.14 
+ -24.14 < Mi <-23.79 
+ -23.79 < Mi <-23.43 
+ -23.43 < Mi <-23.07 
+ -23.07 <Mi<-22.71 
+ -22.71 <Mi<-22.36 
+ -22.36 < Mi <-22.0 



redMaPPer central galaxies 


+ LOWZs 
+ non LOWZs 


N 

O 



‘‘probabilistic-sample'. When using this sample, our calculations 
will account for membership probabilities by summing over Pmcm 
values. 

However, in some cases it is useful to have a deterministic 
sample of LRGs in clusters (e.g., for plotting). In such cases, we se¬ 
lect the A most likely cluster members for each cluster (A is rounded 
to the nearest integer). We will call this the “fixed-sample” We will 
clarify in the text when each of these samples is considered. When 
using the fixed sample, the galaxy with the highest value of Pccn is 
called the redMaPPer central. 


3 RESULTS 

In this section we investigate two questions: (i) How often are cen¬ 
tral galaxies LRGs? and (ii) How often are central galaxies also 
the brightest LRGs in each cluster? For most of our calculations, 
we treat cluster centers and member galaxies probabilistically and 
use the centering probability, Pcen, and the membership probability, 

pmem- 


3.1 Cluster Centroids 

We consider four different cluster centroids: 

(i) RMCG: the redMaPPer central galaxy. For each cluster, this 
is the galaxy with the highest value of Pccn- This is the galaxy that 
redMaPPer selects as the most likely center for each cluster. 

(ii) BMEM: the brightest red-sequence cluster member iden¬ 
tified by redMaPPer, where brightness is given by fcmod- In most 
cases, this is defined with the probabilistic-sample and treated 
probabilistically according to the Pmcm values. In some cases, 
BMEM will be defined using \hs fixed-sample (we clarify in each 
case which definition is used). 


Table 3. Percent of time when two of samples (i) - (iv) overlap. For example, 
in clusters that have at least one LRG, 84.1% of BLRGs are also RMCGs. 
In clusters that have at least one LRG/LOWZ, 94.7% of BLOWZ are also 
BLRGs. Here, BLRG, BLOWZ and BMEM are defined using the “fixed- 
sample”. 



RMCG 

BLRG 

BLOWZ 

BMEM 

RMCG 

100 

84.1 

80.1 

80.5 

BLRG 

- 

100 

94.7 

95.5 

BLOWZ 

- 

- 

100 

96.3 

BMEM 

- 

- 

- 

100 


(iii) BLRG: the brightest DRSPhotLRG in each cluster (if the 
cluster has at least one DRSPhotLRG). This is defined in a similar 
fashion as BMEM. 

(iv) BLOWZ: the brightest PhotLOWZ in each cluster (if the 
cluster has at least 1 PhotLOWZ). This is defined in a similar fash¬ 
ion as BMEM. 

Note that RMCG is defined as the galaxy with the highest 
value of Pcen, but that the actual centering probability of RMCG 
galaxies may be less than 1. In our sample, RMCGs have Pcen val¬ 
ues between 0.27 and ~ 1.00 with a mean value of (Pcen) = 0.87. 
Hence, RMCGs may not always be the true cluster center. Indeed, 
redMaPPer predicts that 13% of RMCGs will not be a true cluster 
center. 

We show a simple comparison between these different classi¬ 
fications in Table]^ We use the “fixed-sample” (see Section [23] for 
the sample definition) to define BLRGs, BLOWZs and BMEMs. 

3.2 Central Occupation Eunctions of LRGs and LOWZ 
Galaxies 

We now compute the probability that a cluster of richness A 
hosts a central galaxy that is also an LRG or a LOWZ galaxy. 
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When expressed as a function of halo mass, this is known as the 
LRG/LOWZ central occupation function, Ncen{M 20 Qb)- To begin 
with, we compute Ween as a function of cluster richness. To mea¬ 
sure Wcen(A), we use the ‘"probabilistic-sample” where the sum of 
Teen for each cluster is equal to 1. 

We use 10 bins in cluster richness from A = 20 to A = 150. 
In each richness bin, we measure the mean value of Pcen for all 
LRG (or LOWZ) members. For any given cluster, the total value 
of Pcen summed over LRG members represents the probability that 
the central galaxy of that cluster is an LRG. In the j-th bin. Ween (A) 
is: 

(Ween(A)), = \ ^ (18) 

where ncius,j(A) is the number of clusters in the j-th richness bin, 
and the subscript i denotes all LRG cluster members in this richness 
bin. 

The results for Ween (A) are shown in the two left panels of Fig- 
ure[^ Erro rs are calculated via bootstrap. For consistency with |Reid| 
|& Sperg^ ( |2009^ , we apply a g-band absolute magnitude cut to our 
LRG sample (we use the “DR8PhotLRG-Mgcut” sample, see Sec- 
tion |2.3| l. For both LRGs and LOWZ galaxies, we find that Wcen(A) 
does not converge to 1 for large values of A. Instead, Wcen(A) flat¬ 
tens and converges to ~ 0.9 at A = 150. 


3.2.1 Non-LRG redMaPPer centrals 


The fact that Ween (A) does not converge to 1 even for large values of 
A can be explained by Figure]^ This figure shows the i-band abso¬ 
lute magnitude distribution of the most likely RMCGs as a function 
of cluster richness. As expected, the absolute magnitudes are cor¬ 
related with richness: brighter galaxies tend to be in richer clusters. 
However, there is significant scatter in this relation. The RMS dis¬ 
persion of Mi as a function of A is (tmi ~ 0.50 for all RMCGs 
and CTMi ~ 0.46 for RMCGs with Pcen > 0.99. The dispersion 
does not depend strongly on A. These values are similar to those 
reported by [Hansen et al.| p005| l for the MaxBCG cluster finder 
Giogj^a(Li) = 0.17, or CTMi = 0.425). In Section|2.2| we show that 


the LRG cuts roughly select galaxies above a fixed absolute mag¬ 
nitude threshold. This is shown by green horizontal lines in Figure 
Because of the trend with Mi and the significant scatter in Mi, a 
sizable fraction (28%) of RMCGs in clusters with A < 40 are sim¬ 
ply too faint to pass the LRG cuts. The scatter in Mi at fixed halo 
mass also means there is a non-negligible fraction of faint central 
galaxies in rich, massive halos that fail to pass the LRG cuts. 

The right hand side of Figure paints a similar picture for 
LOWZ galaxies. The main difference is that the LOWZ selection 
extends to fainter absolute magnitudes so a larger number of RM¬ 
CGs are also PhotLOWZs. There is however, one additional differ¬ 
ence due to the bright r-band cut applied for the LOWZ selection 
(see Equation tn) in Section |2.4^ . This cut excludes about 102 
bright central galaxies (these have a mean r-band magnitude of 
(remod) = 15.8). For some studies, it may be desirable to add these 
bright galaxies back in so as to remove the effects of the LOWZ 
bright cut. Indeed, many of these very bright galaxies are already 
targeted as main sample galaxies or as classical LRgQ Among the 
102 bright RMCGs in our sample that do not pass the LOWZ cut. 


^ This will not be true for parts of the southern galactic cap that were ex¬ 
tended by the BOSS survey. 


73 have spectroscopic redshifts in DRIO. These missing galaxies 
contribute to the low Wcen(A) at high halo mass shown in Figure]^ 


3.2.2 Converting W„„(A) to W„,.(M 2 oo(.) 

We now convert Wcen(A) to Ncm{M 2 Q 0 b) to compare with results 
derived from HOD modeling of counts-in-cylinders and the two 
point correlation function jReid & Spergel||2009| [Parejko et al.| 
|2013| l. In order to convert A to halo mass, we adopt the richness- 
halo mass relation from |Rykoff et al.| ( |2012| l. In Section|^ we per¬ 
form tests to evaluate the impact of the uncertainty in this relation 
on our results. In |Rykoff et al.| ( [2012[ >, the probability that a cluster 
of richness A has mass M 2006 , P(AT 2 ooi)| A), is assumed to follow 
a log-normal distribution with a scatter of aM\x = 0.25. The mean 
of this log-normal relation is given by: 


In 


M 2006 


h^Q 


= 1.72-b 1.08 In ( — 
60 


(19) 


Using this conversion, our lower richness limit of A = 20 cor¬ 
responds to a halo mass of M 2 oob ~ 2.44 x 1C*Mq and a 
richness of A = 100 corresponds to a halo mass of M 2 oob = 
1.39 X lO^^Mg. _ 


Rykoff et al. (2012 > specify P(M 2 ooi)|A) but P(A|M 2 oo 6 ) is 


more useful for our calculations. We derive P(A|M 2 oo 6 ) following 
the method described in the appendix of |Leauthaud et ak|j2010| l 
which assumes that P(A|M 2006 ) is also a log-normal of the form: 


P(A|M2oof,) = 




■ exp 


(In A In Amean {M200b)Y 




( 20 ) 


For this calculation we adopt the halo mass function from | Tinker| 
|et al.|j^08^ . From Equation l |20^ and the halo mass function, we 
find cta|m = 0.231 and 


In Amean (M 200 b) = 0.925 In - 28.26. (21) 

With this relation in hand, we now convert Wcen(A) to Ncen{M 2 Q 0 b). 
Eor this purpose, we assume that Ncen{M 2 oob) follows the func¬ 
tional form: 


{Ncen{M200b)) = 


/a 


1 + erf - 


log^Q(M200b/Afm 


(^log M 


( 22 ) 


This functional form for Ween is similar to that in |Reid & Spergel] 
( |2009^ but with an extra free parameter, /a, which allows Ween to 
converge to values less than unity at large halo masses. 

Given Equation l |20[ > and assuming the [Tinker et al.| p008[ > 
halo mass function ^ at 2 = 0.2, Wcen(A) is computed as : 


^ / P(A|M200b)Ween(M200b)^(M200b)rfM200i, 

“ /P(A|M 200 f,)^(M 2006 )dM 200 f, 

(23) 

We fit our measurements of Ween (A) using Equation l [23[ l with free 
parameters /a and Mmin. The best fit parameters are given in Table 
PI The parameter (TiogAf is fixed to 0.7, which corresponds to the 
value found by [Reid & Spergei]j2009^ . The cluster sample used in 
this paper (A > 20) does not extend to low enough halo masses 
to fully constrain ciog Af- However, the redMaPPer cluster catalog 
does extend to A > 5, which could help constrain ciog m in future 
work. 
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Richness X log (M^oq,, / M, J 


Figure 3. Left panels: Ncen(A) for LRGs (upper panel) and for LOWZ galaxies (lower panel). Solid black lines show our best fit to A^cen(A) when is left 
as a free parameter, and the black dashed lines show the best fit when = Red dash-dotted lines are the results from |Reid & SpergelH2QQ9'l and jParejk^ 
|et al.jpOl^ that have been converted to A^cen(A) using Equation Right panels: A^cen(M2oo6) inferred by fitting Equation W to the measured A^cen(A). 
Solid black lines show our results when /a is a free parameter. The thick black portion of the curve corresponds to the region constrained by the cluster 
data. Black dashed lines show our results when Ja = 1- Red dash-dotted lines correspond to (A^cen) and red dotted lines correspond to (A^sat) from |Reid] 
|& Spergd]|20Q9) (top) and|Parejko et al.|(?013) (bottom). The grey solid line in the lower right hand panel shows A^cen(Af2oo&) we obtain by re-fitting the 
data points from |Parejko et al.H2Q13) . The 68 percent confidence region for this fit is shown using grey shaded regions. Eor both LRGs and LOWZs, we find 
a significant difference between the normalization of Ncen inferred from direct counts in clusters and that from HOD modeling of the two-point correlation 
function (compare black lines to red dash-dot lines). 


Table 4. HOD fits to A^cen(A) for the DR8PhotLRG-Mgcut sample and the PhotLOWZ sample. We fit Ween (A) both fixing fA = l and allowing Ja to be a 
free parameter. The last two columns in this table are derived and discussed in Section]^ 



/a 


X^/dof 

»^CEN [Mpe 3] 

Implied satellite fraction 

LRG 

1.0 (fixed) 

14.19 ±0.01 

15.1/9 = 1.7 

1.8 X 10“® 

47% 

LRG 

0.953 ±0.014 

14.14 ±0.02 

4.87/8 = 0.61 

1.9 X 10-® 

44% 

LOWZ 

1.0 (fixed) 

13.85 ± 0.01 

53.8/9 = 5.9 

3.5 X 10-® 

65% 

LOWZ 

0.947 ± 0.007 

13.69 ±0.03 

15.2/8 = 1.9 

4.2 X 10-® 

58% 
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Figure 4. Left panel: 2 -band absolute magnitude of central galaxies as a function of cluster richness. Red points represent RMCGs that are also DRSPhotLRGs. 
Small dark blue points show RMCGs that fail to pass our photometric LRG selection. Open black diamonds with error bars show the mean magnitude in each 
richness bin with errors that represent the standard deviation of the distribution in each bin. The green shaded region is the approximate 2 -band magnitude 
limit that defines the LRG selection (—23.43 < Mi < 23.07, see Figure^. At A < 40, many RMCGs are simply too faint to pass the LRG cut, and there 
are excluded faint RMCGs at larger richness as well. Right panel: similar plot as in the left panel using the LOWZ selection criteria. Red points represent 
RMCGs that are also PhotLOWZs, and black points represent RMCGs that are not PhotLOWZs. Filled blue triangles are RMCGs that are not PhotLOWZs 
because they are brighter than rcmod = 16. The light-blue horizontal band shows the absolute 2 -band magnitude cut that roughly defines the LOWZ selection 
(-23.07 <Mi< -22.71). 


3.2.3 Comparison with\Reid & Spergel\\2009\ and\Parejko et al.\ 

( |2071t 

In Figurejs] we compare our Ncai{M 20 Qb) with the central occupa¬ 
tion function derived by |Reid & Spergel|p009| l, who used a com¬ 
bination of clustering and the counts-in-cylinders technique to con- 
strain the HOD of LRGs. The halo mass definition assumed inIReid 


& Spergel (2009 i is the same as this worll^ The Reid & Spergel 


1 2009^ results use halos identified by a spherical overdensity algo 


rithm and so will be comparable with the redMaPPer cluster finder. 

For LOWZ, we compare with A'cen derived by|Parejko et al.| 
( |2013r Because |Parejko et al.lpOTs] ) use friend-of-friend halos, 
their results are not directly comparable to ours. We take two mea- 
sures to acco unt for this difference. First, we convert the |Parejto 
et al.H2013| results to our halo mass definitiorp] following [More 


et al.lpoTT). Second, we refit the measurements of abundance and 


clustering of LOWZ galaxies presented by |Parejko et aL|p013) . For 
this, we use the analytical framework for the halo occupation dis¬ 
tribution developed in |van den Bosch et al.H2013) , which accounts 
for the radial dependence of the bias, halo exclusion, and redshift 
space distortion effects on the projected clustering of galaxies (see 


ing data, we use the parameterization as 

Parejko et aL|12013J (see 

also|Miyatake et al.|2013[|More et al.|20 

We adopt a consistent 


® The HOD parameters given in |Reid & Spergel j2009^ assume h = 0.7. 
We have checked that integrating the |Reid & $pergel| j2009) HOD with 
h = 0.7 yields the correct LRG number density. 

® The FOF halo masses roughly conespond to an overdensity of 500 
times the mean matter densi ty jMore et al.|201 1 j . We convert the FOF halo 
masses from |Parejko et ar]j2013t to our halo mass convention assuming 
that M 2 oob /= 1.32. 


halo mass definition of M 2 oob while performing the fits so that the 
HOD modeling results can be directly compared with the results 
from our cluster catalog. The 68 and 95 percent confidence regions 
are shown using grey shaded regions in Figure]^ Our re-fits to the 
data points are in excellent agreement with 
HOD after converting to a single halo mass 
definition, suggesting that differences between friend-of-friend ha¬ 
los and spherical overdensity halos do not have a large impact on 


Parejko et al. 

2013 

Parejko et al. 

(1201311 


iVce 


For both the LRG sample and the LOWZ sample, we find sig¬ 
nificant differences in the amplitude of Ncen derived from HOD 
modeling of LRG/LOWZ clustering and derived directly from the 
redMaPPer cluster catalog. Possible explanations for this discrep¬ 
ancy are presented in Section]^ 


3.3 How often is the Brightest LRG not the Central Galaxy? 


Although many studies assume that the brightest galaxy in a clus¬ 
ter is also th e central galaxy jvan den Bosch et al.|2004[|Weinmann| 
|et al.|[2006{ [Budzynski et al.||2012|>, some recent studies suggest 

that is not always the case (|van den Bosch et al.|2005||Coziol et al.| 

|2009[[Smderson et al.|2009[|Hikage et al.|2013 i. Because redMaP¬ 
Per uses color and position in addition to luminosity to determine 
which galaxy is the central galaxy, we can test this assumption by 
calculating the probability that the brightest LRG/LOWZ galaxy 
in a cluster is not the central galaxy (denoted Pbnc). To compute 
Hbnc. we follow the probabilistic method outlined in Reddick et 
al., in prep. For each cluster, we identify all LRGs brighter than the 
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Figure 5. Probability that there exists a redMaPPer cluster member (or¬ 
ange), an LRG member (black), or a LOWZ member (red), that is brighter 
than the central galaxy as a function of halo mass. Errors are derived via 
bootstrap. Grey points show measurements from Skibba et al.H201 f) . Grey 
lines show predictions from |Skibba et al.] 12011 based on the conditional 
luminosity function (CLF) from Cacciato et al. 2009, and assuming three 
different values for the slope in the satellite CLF (s = 1 dotted cuiwe, 
s = 2 solid curve, s = 3 dashed curve) 


2.8 arcmin 
0.5 M pc 



Figure 6. SDSS image of a cluster in which the redMaPPer central galaxy 
is not the brightest member. Member galaxies are shown with white ticks. 
The richness of this cluster is A = 54.8. The galaxy circled in red in the 
middle is the most likely central, and the one circled in pink (lower left) 
is the brightest member. The centering probability of the central galaxy is 
0.930. The brightest member of this cluster has Pcen = 0.0. The 2 -band 
apparent magnitudes of the central galaxy and the brightest member are 
mi = 16.8 and = 16.3, respectively. 


i-th candidate central. Pbnc for LRG is computed as: 

5 

-Pbnc = E Pcen,z ^ ^ ( 1 Pmem,fc ). (24) 

i=l j k 

The first sum runs over all five candidate centrals, not limited to 
LRGs. The second sum over j is for all LRGs that are brighter than 
central candidate i. The final product over k is for LRG members 
brighter than galaxy j to ensure that only the BLRG is considered. 

Figure shows the mean probability that the BLRG or 
BLOWZ galaxy in a cluster is not the central galaxy. We find that 
Pbnc varies between 0.2 and 0.3 and mildly increases with halo 
mass. Figure]^ shows an example of a cluster in which the central 
galaxy is not the brightest member. The most likely central galaxy 
of this cluster has Pcen = 0.930, and the Pcen value for the bright¬ 
est galaxy (in the, fixed-sample) is 0. The central is not an LRG but 
the brightest member is an LRG. In this cluster, the BLRG is less 
likely to be the central because it is offset from most other cluster 
members. 

We compare our results with [Skibba et al.| (|2011) in Fig¬ 
ure [Skibba et al.| ( |201 1[ | compute Pbnc by comparing the ve¬ 
locity and positional offsets of the brightest cluster galaxy from 
the other cluster members to the expected offsets computed using 


mock group catalogs. The Skibba et al.| ( |20ir| l results are derived 
from the [Yan g et al.jpOO? i group catalog. The details of the clus¬ 
ter finding algorithm are accounted for in thejSkibba et al.|(2011^ 


analysis by running thel^n^e^l. 

12007 1 group finder on mock 

catalogs. Note that thejSkibba et a 

L 

( 2011 ) results span a lower 


redshift range ( 0.01 < z < 0 . 20 ), and larger halo mass range 
(12 < log(Mhaio/M 0 ) < 15) than our results. We convert the 
[Skibba et al.| ( |20lT| l results to our assumed value of h = 0.7. We 
do not, however, account for the difference in halo mass definition, 
Misob in [Skibba et al.[ ( [T01 l[ l versus M 2 oob here, as the correction 
is of order 3%. 

Our values of Pbnc are significantly lower than those found by 


Skibba et al.ji 2011 1 . At a halo mass of log(Mi 8 of)/AT 0 ) = 14.29, 


Skibba et al.H20lT) find Pbnc = 0.43, whereas we find Pbnc = 


0.23 at log{M 2 QQb/M q) = 14.5. As demonstrated in Reddick et 
al., in prep, part of the difference may be explained by a tendency 
to over-subtract the sky around bright, extended galaxies in SDSS 
data releases prior to DR 8 . While we rely on DR 8 data,[Skibbaetal.| 
( [2011^ use SDSS DR4 photometry. Improvements in the sky sub¬ 
traction in DR 8 typically increased the measured luminosities of 
central galaxies more than those of satellites ( [Bemardi et al.|20T3) >, 
thus decreasing the probability that a satellite is brighter than the 
central galaxy. However, the sky subtraction is most relevant at red- 
shifts below the range we study here, and may only explain part of 
the discrepancy. We refer to Reddick et al., in prep for further de¬ 
tails on this question. 

In conclusion, we find that within the redMaPPer cluster cat¬ 
alog, the brightest LRG/LOWZ galaxy is also the central galaxy in 
only ~70% of clusters. This will impact redshift space distortions 
(RSD) studies that construct “halo catalogs” under the assumption 
that the brightest LRG is always the central galaxy ( [Hikage & Ya- [ 
[mamoto|2013^ . This also suggests that naive cluster-finding algo¬ 
rithms that assume the central galaxy is the brightest galaxy in a 
cluster may have a 20 ~ 30% mis-centering fraction ivan den| 


[Bosch et al.|2005][Coziol et al.|200^[Sanderson et al.|2009|. 


3.4 Offset Distributions between the Central Galaxy and the 
Brightest LRG 

We now investigate how the radial distributions of non-central 
BLRGs compare with the radial distributions of red cluster mem¬ 
bers. We compute the distribution of projected offsets between the 
brightest and central galaxies for clusters in which the BLRG is not 
the central galaxy. Specifically, we compute Doff/P a, where Doff 
represents the 2D projected distance between the central galaxy and 
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Offset / 

Figure 7. Distribution of _Doff/R a for the brightest LRGs, when the RMCG 
and BLRG are different. Errors are calculated by bootstrap. The red solid 
curve shows the best fit assuming a Rayleigh distribution. Black crosses 
represent the distribution of all redMaPPer cluster members. Black plus 
signs represent the distribution of all redMaPPer cluster members in clus¬ 
ters in which the BLRGs are not the central galaxies. The fact that the 
black crosses and black plus signs trace similar distributions suggests that 
clusters with non-central BLRGs are not dominated by projection effects 
and/or mergers. Grey curves represent the expected distribution of satel¬ 
lites distributed according to NEW profiles with different halo masses 
at redshift z = 0.25. Grey solid, dashed, and dash-dotted curves are 
log(AT 200 i)) = 14.0,14.5,15.0, respectively (and the halo concentrations 
are c = 7.2, 6.5, 5.3). The distribution of D^ff/Rx for non central BLRGs 
is significantly different (truncated at the outskirts) compared to the distri¬ 
bution for all cluster members. 


bution of radial distances to all cluster members. The differences 
between these distributions may indicate a larger impact of dynam¬ 
ical friction on more massive sub-halos as suggested hy |Wu et al.| 
( |2013^ . Figure 7 also shows that the distribution of radial offsets for 
member galaxies in clusters is the same whether or not the BLRG 
is the central galaxy. This suggests that clusters with non-central 
BLRGs are not dominated hy odd (e.g. merging) clusters. 

Previous work has sometimes modeled the distribution of off¬ 
sets using a Rayleigh distribution function instead of an NFW pro¬ 
file ( [Johnston et al.|2007l ): 


P{d) = Az exp 
cr^ 


1 

2 


2 


(26) 


where d = Dofi/Rx- As shown in Figure|7] a Rayleigh distribution 
function with a — 0.39 (our best fit value) does not provide a 
satisfactory description of this distribution, although it provides a 
better fit than an NFW distribution. 


4 DISCUSSION ON THE ORIGIN OF THE Wcen 
DISCREPANCY 

The above results show that (i) the fraction of redMaPPer clus¬ 
ters with a central LRG (Acen(A)) does not converge to unity for 
rich clusters, (ii) the brightest LRG in a cluster is not the central 
galaxy in 20 ~ 30% of cases, and (iii) the offset distributions of 
non-central BLRGs is significantly different than that of all cluster 
members. Below, we put forth possible explanations for why our 
measured Nan is lower than that inferred from clustering. 


the BLRG/BLOWZ. We normalize Dog by the cluster cutoff radius, 
Rx, to account for size variations among clusters. Dott is computed 
as follows. For each cluster, we consider all 5 central galaxy can¬ 
didates. Let us consider the i-th central candidate with a centering 
probability of Pan,i- We identify all LRG cluster members brighter 
than cenhal galaxy i. The probability of observing an offset Doff,*i 
between central galaxy i and the brightest LRG j in the cluster is: 

P{Dog,ij) = Pcnn,iPmem,j Pmem,k)j (25) 

k 

which is simply the i, y-th term in the sum making Pbnc (see Equa¬ 
tion]^. 

Figure]^ is constructed by finding all configurations in which 
LRG j is brighter than central galaxy i and computing Dog,ij and 
P{Doif,ij)- We then bin the Dns/Rx values and sum the P{Dns,ij) 
values in each bin. We only show the distribution of offsets for 
BLRGs, because the distributions are similar for BLOWZ galaxies 
and BMEM galaxies. The mean offset between the RMCG and the 
brightest LRG (when these galaxies differ) is 0.4i?A- 

For reference, we show the expected distribution of Dos/Rx 
assuming an NFW profile (grey lines). This is computed by first 
integrating the three dimensional density profile along the line-of- 
sight to obtain the projected surface mass density. The curves are 
then computed by multiplying the area of a thin annulus at each ra¬ 
dius to the density at each radius. Each grey line shows a different 
halo mass and halo concentration. Eigure shows that the distri¬ 
bution of 2D projected halo-centric distances to BLRGs is signif¬ 
icantly different (more truncated at the outskirts) than the distri- 


4.1 Possible (but unlikely) explanations 


One reason that our values for Nan are lower than those found in 
previous studies might be that the redMaPPer algorithm does not 
correctly identify central galaxies. We have conducted a visual in¬ 
spection of 272 redMaPPer clusters in our redshift range. We cre¬ 
ated color images for these clusters using the Stripe 82 coadds jAn-j 
jnis et al.|2014^ which are ~ 2 mags deeper than the SDSS single 
epoch imaging. Eive co-authors inspected these clusters to select 
“visual” central galaxies. Our visual inspection did not reveal any 
obvious failure modes in the redMaPPer cenhal galaxy selection 
and we found good agreement between the “visual” central galax¬ 
ies and the redMaPPer galaxy with the highest value of Pan (95% 
of the time). From these tests alone we cannot rule out the possibil¬ 
ity that the differences that we find for Nan are due to an incorrect 
central galaxy assignment by redMaPPer. However, our visual in¬ 
spections did not reveal any obvious issues. In Paper II we perform 
tests of the redMaPPer centering probabilities using weak gravita¬ 
tional lensing and cross-correlations. Early results from this work 
suggest that the redMaPPer centering scheme outperforms other 
centering choices, such as the bright cluster member. 

Based on the above tests, we think that mis-centering in the 
redMaPPer cluster catalog is unlikely to explain the Nan discrep¬ 
ancy. Another possibility is that part of the difference may be 
caused by uncertai nties in th ejRykof f et al.| ( |2^012| l mass-richness 
relation. Indeed, the jRykoff et ah j20I2^ m^-richness relation was 
derived for an older version of the redMaPPer cluster catalog. To 
determine how large of an impact this might have, we use an up¬ 
dated calibration from Reddick et al. in prep, which is based on 
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the same redMaPPer catalog as the one used herj^ We re-derive 
our results using the updated mass-richness relation but as shown 
in Figure this only slightly alters our results. Finally, we also 
vary aM\\ between aM\x ~ 0.22 and aM\x ~ 0.28 but this does 
not impact our conclusions. Based on these tests, we conclude that 
uncertainties in the |Rykoff et al.| ( [2012^ mass-richness relation are 
unlikely to explain the discrepancy. 

Another possibility is that there is a systematic difference be¬ 
tween our photometric LRG selection and the actual spectroscopic 
samples. However, we obtain similar values for NccnW if we se¬ 
lect LRGs based on DR7 photometry, or based on the original LRG 
targeting flag (TARGET_GALAXY_RED, see Appendix A). As dis¬ 
cussed in the Appendix A, differences between our photometric 
LRG selection and TARGE T_GALAXY_RED are consistent with 
scatter in the photometry between reductions. Hence, differences 
between the LRG samples from |Reid & S pergel H2009 ^ and lPare-] 
[jko et al.| ( |2013| l and our photometric LRG selection are unlikely to 
be the source of the Ncen discrepancy. 

Another possibility may be related to issues in our conversion 
from Acen(A) to Acen(M 2 ooi)). For example, for consistency with 
|Reid & Spergel|p009^ , we fixed the parameter aiog m to 0.7. How¬ 
ever, we have also tested fits with criogM = 0.5 and ciogM = 0.9. 
The shape of NcsniM 2 oob) is modfied as expected but the overall 
discrepancy persists. The conversion of Wccn(A) to Ncen{M 2 oob) 
also depends on the assumption that P{M 2 oob\^) follows a log¬ 
normal distribution. While this is a very common assumption, it is 
by no means a well tested one. It would be interested to test how 
this assumption would affect our results, but we defer this aspect 
to future work. Finally, this conversion may also fail if redMaPPer 
clusters do not represent a perfect recovery of dark matter halos. If 
the cluster catalog contains purity and/or completeness issues, then 
it will not be safe to assume that one can derive Ncen{M 2 oob) from 
/Vcen(A). Completeness, however, does not appear to be an issue for 
the redMaPPer cluster catalog, beyond the expected scatter in the 
richness-mass relation, which we have already accounted for in the 
density matching procedure. Purity however, is more of a concern. 
Preliminary X-ray follow-up of a complete sample of redMaPPer 
clusters suggests that the fraction of redmapper clusters that are 
projections is no larger than 7%. A 7% change in density corre¬ 
sponds to a ~2% change in mass which corresponds to a ~2% 
change in the scaling relation ( [Weinberg et al.|[2013^ and should 
not dramatically alter our conclusions. 


4.2 Implication of our Results in Terms of Central and 
Satellite Number Densities 

In order to compare our results to previous work based on clus¬ 
tering, we compute the number density of central and satellite 
galaxies implied by our measurements of Wcen(A). To compute 
the central number density, we integrate the product of the best- 
fit Ncen{M 2 oob) and ■^(M 2 oob) over the full halo mass range. 
For this calculation, we are forced to extrapolate our model for 
Ncen{M 2 oob) to lower halo masses that are not sampled by the 
redMaPPer galaxy clusters used here. There is a potentially large 
systematic uncertainty associated with this extrapolation. On the 
other hand, the erf model in Equation (|22[ is the standard for 


Because Reddick et al. in prep assume a different cosmology than 
we do, we recalibrate their mass-richness relation for this test. This re¬ 
calibration yields the same cluster abundances as a function of A as Reddick 
et al. in prep. 


HOD models. Therefore, any discrepancies between our work and 
clustering studies that can be traced back to the extrapolation of 
Ncen{M20Qb), would Correspond to a failure in the now standard 
HOD model. Thus, we simply caution that our conclusions are de¬ 
pendent on the validity of the erf extrapolation. Also note that the 
goal of this exercise is simply to gain a sense of the differences 
between our results and previous work based on clustering. 

Our HOD implies a number density of central galaxies of 
1.9 X 10“® Mpc“® for LRGs and a number density of 4.2 x 
10“® Mpc“^ for LOWZ galaxies (see Table 0- Assuming a to¬ 
tal number density of titot = 3.4 x 10“® Mpc“^ for LRGs and 
titot = 1.0 X 10“'* Mpc“^ for LOWZ galaxies, we can infer the 
number density of satellite galaxies (tisat = utot — ncEN) and the 
implied satellite fraction (see Table|^. For LRGs we find that our 
model implies a satellite number density of 1.5 x 10“® Mpc“® 
which corresponds to an LRG satellite fraction of 44%. Likewise, 
for LOWZ, we find that our model implies a satellite number den¬ 
sity of 5.8 X 10“® Mpc“^ which corresponds to a satellite fraction 
of 58%. These satellite fractions are much higher those ofjR eid &| 
|Spergel|p009| l and jPareJko et al.|p013| (, which are 6.4% and 12%, 
respectively. 

The fact that our HOD fits to Nain{M 2 oob) imply much larger 
satellite fractions than previous works is not surprising because our 
Ncen{M 20 Qb) is lower than those of Reid & Spergel i2009i and 
jPareJko et al.H20i^ . Hence, the satellite fraction must increase in 
order to preserve the total LRG number density. The satellite frac¬ 
tions derived in this section depend on multiple assumptions and 
on extrapolating below our halo mass limit. Nonetheless, the satel¬ 
lite fractions inferred from our fits are uncomfortably high — lead¬ 
ing us to explore effects that might invalidate our measurements of 
Ncen{M 20 Qb)- Since our measurements of A/;en(A) are reliable (see 
§4.1| l, we focus on effects that might influence the conversion from 

Wcen(A) to Wcen(M2006). 


4.3 How would our Results be Affected by a Strong 

Correlation Between Richness and Central Luminosity? 

Our results rely on converting A to M 2 oob via Equation jl9| l, which 
assumes that A only depends on M 2 oob- Here, we consider the pos¬ 
sibility that A also correlates with Mp®, the absolute magnitude of 
the central galaxy. 

In this case, the conversion in Equation l |19[ ( will depend on 
P(A|M 2 ooi., Mp^) instead of just P(A|M 2 oofj). To determine how 
this correlation might impact our results, we evaluate Ncm{M 2 oob) 
for two mock samples in which halo mass and central galaxy lumi¬ 
nosity are maximally correlated and anti-correlated, respectively, at 
fixed richness. We design these samples in the following way: 

(i) We build a mock catalog that contains halos drawn from the 
[Tinker et al.| ( |200^ mass function. Each halo is assigned a richness 
based on Equation ( [20[ l. 

(ii) We assign an i-band central galaxy absolute magnitude to 
mock halos. Mock Mp® values are randomly drawn from the real 
redMaPPer catalog after matching mock clusters and real clusters 
by richness. This procedure ensures that our mock catalog has the 
same relationship between A and as the redMaPPer catalog. 

(iii) We then bin the clusters by A. Each narrow richness bin 
contains a range of halo masses and central galaxy magnitudes. 
In the initial catalog, there is no correlation between Mp^ and 
M 2 oob- Within each bin, we will reassign the Mp^ values in order 
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to introduce a correlation between halo mass and central galaxy 
magnitude. 

(iv) In our first model (MODI), we reassign all central galaxies 

such that heavier halos contain brighter centrals within each rich¬ 
ness bin (i.e., M 2001 . and are maximally correlated at fixed 

A). In this model, clusters with lower values of A have brighter cen¬ 
tral galaxies at fixed halo mass (i.e., A and are anti-correlated 
at fixed halo mass). 

(v) In our second model (MOD2), we reassign all centrals galax¬ 

ies such that heavier halos contain fainter centrals at fixed richness 
(i.e., M 2 oob and are maximally anti-correlated at fixed A). 

In this model, clusters with lower values of A have fainter central 
galaxies at fixed halo mass (i.e., A and Mp^ are correlated at fixed 
halo mass). 


We compute Ncen{M 2 oob) for MODI and MOD2 and the re¬ 
sults are shown in Figure]^ Interestingly, we find that correlations 
between halo mass, richness, and central galaxy luminosity may 
have a large impact on our inferences about A^cen(M 2 oo 6 ). If the 
magnitude of the central galaxy is maximally correlated with halo 
mass at fixed richness, then our results would agree with those of 
|Reid & Spergell ( |2009| l. For LOWZ galaxies, however, there is a 
conflicting effect due to the bright cut given by Equation l |15^ . This 
cut removes bright galaxies, which, in MODI, are preferentially lo¬ 
cated in the most massive halos. Therefore, Ncen{M 20 Qb) decreases 
at high halo masses for MODI. By removing the bright cut, we ob¬ 
tain similar results for the LOWZ and classical LRG samples. How¬ 
ever, a fair comparison to |Parejko et al.| ( [2013) requires the bright 
cut for LOWZ selection. 

A strong anti-correlation between central luminosity and clus¬ 
ter richness (MODI) at fixed halo mass would bring our results 
into good agreement with Reid & Spergel 1 2009 1 and Ncen{M 20 Qb) 
would converge to unity at large halo masses. It is possible that 
merging might lead to such an anti-correlation. When a satellite 
merges with the central galaxy, the richness of the cluster decreases 
while the magnitude of the central galaxy increases. In this sce¬ 
nario, halo mass estimators based on combinations of Mp^ and A 
would result in a smaller scatter than halo mass estimators based 
on just A. In fact, [Reyes et al.| ( (2008| l present evidence this may be 
the case using the MaxBCG cluster catalog. 

Figure demonstrates that correlations between central lu¬ 
minosity and cluster richness would clearly impact our inference 
about Ncea{M 2 oob)- Figure |9| also demonstrates that such correla¬ 
tions would also impact our inference of Pbnc- To match the jReidj 
|& Sperg^ ( |2009[ result, we must assume that central galaxy lumi¬ 
nosity and cluster richness are strongly anti-correlated. While such 
a strong anti-correlation seems unlikely, preliminary results from 
Reddick et al. in prep based on an analysis of the Conditional Lu¬ 
minosity Function (CLF) of redMaPPer cluster members also sug¬ 
gest a strong anti-correlation between richness and central galaxy 
luminosity. 

In conclusion, a strong anti-correlation between central lumi¬ 
nosity and cluster richness at fixed halo mass is required in order to 
reconcile our results with those based on clustering studies. How¬ 
ever, further investigation is needed to confirm or repudiate this 
hypothesis. An interesting direction for future work will be follow¬ 
up on our study by performing a joint CLF analysis of LRG clus¬ 
tering, galaxy-galaxy lensing, and direct counts in clusters (|Tinker| 
jet al.|2012| l. 
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Figure 9. Impact of MODI and MOD2 on our inference of Pbnc- Black 
plus signs, blue and magenta diamonds represent PfeNC values when there 
is no correlation between central luminosity and cluster richness, and cases 
of MODI and MOD2, respectively. Correlations between central luminosity 
and cluster richness clearly have a large impact on Pbnc - The MODI model 
brings our results into closer agreement with| Reid & Spergel| [2009[ but 
leads to a larger tension with|Skibba et al.H2011[. 


4.4 How Robust is Acen Inferred from Clustering and CIC 
Studies? 

Finally, we note that the values of Acen from studies similar to|Reid| 
|& Sperg^j2009^ that use observables such as counts-in-cylinders 
or clustering could be dominated by theoretical prejudices. For ex¬ 
ample, at the low mass end, HOD models often make the standard 
assumption that a halo cannot host a satellite galaxy if it does not 
also host a central galaxy from the same sample. Under this as¬ 
sumption, most LRGs in low mass halos (which contain on average 
only one LRG or less) will be centrals. But given that observa¬ 
tions of counts-in-cylinders probe Aot they should not be able dis¬ 
tinguish between this standard scenario and a scenario in which a 
fraction of LRGs in low mass halos are in fact satellites. 

The clustering signal may also be relatively insensitive to such 
effects. In these small halos there are no central-satellite or satellite- 
satellite pairs of galaxies. Therefore, the difference in the clustering 
between the two scenarios will arise only due to galaxy pairs in 2 
different halos (the 2-halo term), which also may not be heavily af¬ 
fected. In particular, if non-central LRGs in low mass halos follow 
a radial distribution peaked towards the center as shown in Figure^ 
the differences in the clustering signal may be small. 

Similar arguments hold at the high mass end. HOD models 
often assume that Acen approaches unity at large halo masses. At 
these mass scales, the one halo term in the clustering signal is dom¬ 
inated by the satellite-satellite term (which is not sensitive to Acm), 
except possibly at the smallest scales where the observational un¬ 
certainties are significantly larger jMasjedi et al.|2006l >. As argued 
before, counts-in-cylinders will not be able to make a distinction 
between centrals and satellites either. 

The effect of the above theoretical prejudices on the Acen 
constraints obtained from studies relying on counts-in-cylinders or 
clustering can be studied in detail using mock galaxy catalogs. Such 
an investigation, however, is beyond the scope of the current paper, 
and is left for future work. 
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Figure 8. Tests to explain the origin of the A^cen(Af2006) discrepancy. Our original results from Figurej^are shown with black/grey solid curves. First, we 
test the impact of the mass-richness relation. The dotted grey lines shows A^cen(Af2006) inferred using an updated mass-richness relation with a slope of 1.2, 
instead of 1.08, in Equation (19). The updated mass-richness relation only has a minor impact on our results. Second, we test if correlations between halo 
mass, richness, and magnitude of the central galaxy may impact our results. In both panels, the blue data points show MODI, in which halo mass and central 
galaxy luminosity are correlated at fixed richness (or brighter centrals live in less rich clusters at fixed halo mass). The magenta points show MOD2, in which 
halo mass and central galaxy luminosity are anti-correlated at fixed halo mass. In the case of MODI, the bright LOWZ cut removes the brightest centrals 
which are in the most massive halos causing Ncen{M 2 oob) to decrease in higher mass halos. Grey crosses show MODI when the bright cut is not applied - 
demonstrating that the down-turn in indeed caused by the bright cut. The red cuiwes in the left and right panels show clustering results from|Reid & Spergel| 
j2009) and|Parejko et al.H2013), respectively. 


5 SUMMARY AND CONCLUSIONS 

In this paper, we investigate the connection between central galax¬ 
ies in clusters and two classes of red luminous spectroscopic galax¬ 
ies: “classical LRGs”, which were targeted as part of the SDSS-I 
and SDSS-II programs ( [Eisenstein et al.|2001^, and LOWZ g alax- 
ies from the SDSS-III BOSS program ( [Dawson et al.|2013^ . The 
aim of this paper is to address two different questions: (i) what 
fraction of cluster central galaxies are classified as either an LRG 
or a LOWZ galaxy, and (ii) when a central cluster galaxy is either 
a LRG or a LOWZ galaxy, how often is the central galaxy also the 
brightest LRG in the cluster. 

To tackle these two questions, we use the state-of-the art 
redMaPPer cluster catalog in the redshift range 0.16 < z < 0.33 
and in the richness range A > 20, corresponding to a halo masses 
above 2.44 x 10^'^Mq. To avoid complications due to fiber colli¬ 
sions, we construct a set of photometric LRG/LOWZ cluster mem¬ 
bers. Extensive tests show that these photometric samples are in 
excellent agreement with the original LRG and LOWZ target cata¬ 
logs (> 95%, see Appendix). 

A key feature of the redMaPPer cluster catalog is that it de¬ 
fines cluster members and central galaxies probabilistically. Using 
the centering probabilities given by redMaPPer, we derive the cen¬ 
tral occupation of LRGs and LOWZs as a function of richness, 
Wcen(A), and halo mass, Wcen(M 2 oo 6 ). We find the striking result 
that the central occupation function for both classical LRGs and 
LOWZ galaxies does not converge to unity, even for large clusters. 
Instead, Ween (A) converges to ~ 0.95. This can be explained by the 
large scatter in central galaxy luminosity at fixed richness and the 
LRG/LOWZ selection cuts, which both impose absolute magnitude 
limits that exclude central galaxies in the redMaPPer catalog. 


When we naively convert Wcen(A) into Na:n{M 2 oob), assum¬ 
ing that the halo mass only depends on cluster richness, the central 
occupation still only converges to 0.95 for massive halos. This is 
in conflict with the assumptions behind most HOD-type studies of 
the two-point correlation function. 

While inaccuracies in the redMaPPer centroids could cause 
this discrepancy, we will show in Paper II that the redMaPPer cen¬ 
tering probabilities are reliable. Instead, we show that if the magni¬ 
tude of the central galaxy is maximally correlated with halo mass 
at fixed richness (i.e., A and Mp^ are anti-correlated at fixed halo 
mass), then our results for “classical” LRGs would agree with those 
of |Reid & Spergel]p009| ). This maximal anti-correlation does not 
alleviate the tension between our results for LOWZ galaxies and 
those of |Parejko et al.| ( |2013| l. This is because the LOWZ selec¬ 
tion includes a bright cut, and, by correlating galaxy brightness and 
halo mass, we systematically lower the fraction of central LOWZ 
galaxies in massive halos, thus decreasing Wcen(M 2 oof)) for mas¬ 
sive clusters. 

Our results may indicate a strong correlation between cen¬ 
tral galaxy luminosity and halo mass at fixed cluster richness. This 
is equivalent to a strong anti-correlation between central galaxy 
brightness and richness at fixed halo mass. It is possible that merg¬ 
ing might lead to this anti-correlation: mergers between satellite 
galaxies and the central galaxy will increase the central galaxy lu¬ 
minosity and decrease the cluster richness. In fact, examples of the 
extreme cases of such a correlation have been well studied and 
are known as “fossil groups” (e.g, |Jones et d^|2003t |D’Onghi^ 
|et al.|2003]|Zentner et al.|2005^ . However, further investigation will 
be necessary in order to confirm if this correlation may be strong 
enough to fully explain our results. Another possibility is that the 
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HOD models from counts-in-cylinders or clustering may be sub¬ 
ject to theoretical prejudices in this high halo mass regime. It is 
possible that the differences that we observe for the central occupa¬ 
tion function may be caused by a combination of both effects. An 
interesting direction for future work will be to try to disentangle 
these effects by performing a joint CLF analysis of LRG cluster¬ 
ing, galaxy-galaxy lensing, and direct counts in clusters as in [Tin- 1 
|ker et dL\\\2Q\2) . 

We also investigate how often the brightest cluster member is 
not the central galaxy, Pbnc- Using both the cluster member prob¬ 
abilities and the central galaxy probahilities from redMaPPer, we 
find that in 20 ~ 30% of clusters the brightest galaxy is not the cen¬ 
tral galaxy. Our estimate is in good agreement with a CLF analysis 
by Reddick et al. in prep, but our measurement of Pbnc is nearly a 
factor of two lower than that of Skibba et al. ]M- The large dis¬ 
crepancy with |Skibba et al.| ( |201 1 1 may be in part due to improve¬ 
ments in the SDSS sky subtraction between DR4 (used by|Skibba| 
|et al.|20lT| l and DR8 (used here). Reddick et al. in prep demonstrate 
that improved sky subtraction around luminous, extended galaxies 
tends to increase the difference in luminosity between cenhal and 
satellite galaxies, thus substantially decreasing Pbnc- However, fur¬ 
ther testing is required to determine if the differences in photometry 
apply to galaxies in the redshift range studied here. We have also 
quantified the radial distribution of bright non-central LRGs and 
found that they follow a substantially different distribution (tmn- 
cated at the outskirts) compared to the distribution for all cluster 
members. 

Because we find a significant fraction of clusters in which 
the BLRG is not the central galaxy, using the brightest galaxy 
as a tracer for the cluster center will lead to errors due to “mis- 
centering,” particularly in RSD studies based on LRG “halo cat- 
alogs” (e.g.lReid et al.|20ldl|Hikage et al.|2013) . The analysis in 
[Hikage et al.H2013| suggests a mis-centering fraction for BLRGs of 
30 ~ 40%, in rough agreement with the 20 ~ 30% mis-centering 
we find above. [Hikage et aL] ( |2013| ) demonstrate that mis-centering 
of this magnitude may significantly effect cosmological parame¬ 
ter estimation in future surveys, such as eBOSS, PFS, and DESI. 
Future work will investigate how well the BLRG satellite popula¬ 
tion can be characterized using a joint analysis of BLRG clustering 
and information from cluster catalogues, such as redMaPPer. By 
characterizing satellite BLRGs, the effects of mis-centering may 
be mitigated. 

The results presented here rely on the assumption that the 
probability values from redMaPPer, particularly those for the clus¬ 
ter center, are accurate. A comparison of X-ray to optical centers 
in galaxy clusters with available data suggests that the redMaPPer 
centering probabilities are indeed accurate ( |Rozo & Rykoff]2014| >. 
However, the sample sizes for these comparisons are still small, 
and rely heavily on X-ray selected sub-samples of clusters. More 
detailed investigation into the redMaPPer centering probabilities is 
warranted. Following the methods of |George et al.|j2012^ an d |HlF| 
|age et al.]j2013^ , in Paper II, we will perform tests of the redMaP¬ 
Per centering probabilities by using a combination of weak grav¬ 
itational lensing and the projected correlation between the LRG- 
inferred halos and a fainter photometric sample of galaxies. Pre¬ 
liminary results from Paper II indicate that the redMaPPer centroids 
are better tracers of halo centers than the brightest cluster members. 
This lends credence to both the redMaPPer centering probabilities, 
and the results presented in this work. 

In conclusion, LRGs are considered among the best under¬ 
stood samples of galaxies and are expected to occupy halos in a 
relatively simple way. However, in this paper we have shown that 


even these relatively simple galaxies still harbor suprises and that 
the connection between LRGs and dark matter halos may not be 
straightforward. 
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APPENDIX A: OUR PHOTOMETRIC SAMPLE AND SDSS 
TARGET SELECTION 

Our LRG sample selection is described in Section [2^ Here, we 
compare this sample to the actual LRG sample targeted for spec¬ 
troscopy in SDSS ( [Eisenstein et al.||2001) . Although the cuts de¬ 
scribed above are the same as those for the targeted sample, the 
SDSS photometry changed slightly between DR7 and DR8. As we 
show below, this slightly alters our LRG sample when compared 
to the targeted sample. To do this comparison, we introduce the 
following new samples: 

(i) DR7PhotLRG: Photometrically selected LRGs with DR7 
photometry 

(ii) TargLRG: Original set of LRGs that are targeted by SDSS 
(hi) TargLOWZ: Original set of LOWZs that are targeted by 

SDSS 


Al DRSPhotLRG, DRTPhotLRG and TargLRG 

The SDSS LRGs selected for spectroscopy are flagged 
TARGET_GALAXY_RED, and we queried DR7 to get all 
such galaxies. We call these SDSS LRGs “TargLRGs”. These 
galaxies are selected using the cuts described in Section [23] but 
using earlier reductions of the SDSS data. 

We note here that there is an implicit bright cut in TargLRG, 
due to the bright cut in the SDSS main sample, which is not part 
of our photometric LRG selection. The bright limit is included 
because very bright objects may cause contamination of spectra 
or saturation of the CCDs. Galaxies with apparent magnitudes 
brighter than 15 in g- or r-band or 14.5 in i-band, or Petrosian 
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magnitudes brighter than rpcuo < 15 are excluded from SDSS main 
sample, and, therefore from TargLRG. However, since every galaxy 
in our redMaPPer cluster catalog is fainter than this limit, the bright 
limit for classical LRGs is not relevant to our work. 

We compare TargLRG to our sample of photometrically se¬ 
lected LRGs, DRSPhotLRG. Since the DR8 sky coverage is wider 
than that of DR7, we do the comparison in a region where DR8 and 
DR7 overlap, i.e., the north Galactic cap. However, we find that 
the TargLRG sample contains holes in this region, due to missing 
runs. Therefore, we restrict our comparison to smaller region in the 
north Galactic cap: 120 < ra < 250 and 20 < dec < 50. In this 
region, there are 3,549 DR8PhotLRGs and 3,542 TargLRGs; 3,191 
galaxies are in both samples, which is 89.9% of the DR8PhotLRG 
sample and 90.1% of TargLRG sample. 

We perform the same comparison, replacing DR8PhotLRGs 
with DR7PhotLRGs. These samples are selected identically, but 
the latter uses DR7 photometry. Since the LRG target catalog was 
made before DR8, and there are small changes in the photometry 
between DR8 and earlier data releases, it is worthwhile to compare 
TargLRGs with DR7PhotLRGs. Here we again take the region lim¬ 
ited by 120 < ra < 250 and 20 < dec < 50 in the north Galactic 
cap so as not to be affected by missing runs in TargLRG. In this re¬ 
gion, there are 3532 DR7PhotLRGs and 3542 TargLRGs. The sam¬ 
ples overlap for 3323 galaxies, which is 97.9% of DR7PhotLRG 
and 97.7% of TargLRG. 

In conclusion, when we use DR7 photometry, we find a 6% 
difference between the original LRG target catalog and our photo¬ 
metric selection. When we use DR8 photometry, the difference in¬ 
creases to 12%. This demonstrates, that while similar, our classical 
LRG sample is not identical to the SDSS spectroscopic LRG sam¬ 
ple and care must be taken when comparing the two. The source of 
these differences is likely scatter in the photometry between reduc¬ 
tions. Directly comparing DR7 and DR8 photometry for the same 
galaxies, we find the RMS scatter in the magnitudes and colors 
is ~ 0.07 mags, with no systematic offsets. Photometric scatter 
will affect faint LRGs more than bright LRGs. Reassuringly, when 
comparing TargLRG, DR7PhotLRG, and DRSPhotLRG, galaxies 
in one sample but not the others are on average half a magnitude 
fainter than typical LRGs. Photometric scatter mainly affects the 
LRG sample near the color-magnitude cuts for LRG selection, but 
the overall number of LRGs selected in each sample is nearly un¬ 
changed. Therefore, while the samples of LRGs cannot be directly 
compared, they are equivalent for our analysis. 


A2 BOSS LOWZ target selection 

As with classical LRGs, we compare our original PhotLOWZs 
to the LOWZ target catalog. For the LOWZ target catalog, we 
use “bosstarget-lrg-main007-collate.fits”P^ This catalog contains 
LOWZ galaxies, as well as CMASS-selected galaxies at higher red- 
shifts (0.4 < 2 < 0.7). We select LOWZ galaxies with the flag 

(bossjargetl AND 2°) NE 0. 

We call the LOWZs in the target catalog “TargLOWZs”. Before 
comparing TargLOWZs and PhotLOWZs, we need to remove re¬ 
gions with problems in the target catalog. Due to a bug in the star- 
galaxy separation used to create the target catalog, there are some 


http://www.sdss3.org/drl0/algorithms/boss_target_selection.php 


regions we need to exclude. To remove undesired regions, we re¬ 
move objects from TargLOWZ with 

TILE < 10324. 

Then we simply take the regions that are not affected by this bug. 
They are: 

• 0 < RA < 50, 2 < DEC < 40 

• 140 < RA < 270, 45 < DEC < 70 

• 100 < RA < 220, 6 < DEC < 28 

• 310 < RA < 360, -15 < DEC < 40 

There are 12,778 PhotLOWZs in these regions, 12,294 of 
which (96.2%) are also selected as TargLOWZ. As with classi¬ 
cal LRGs, the photometric LOWZ selection closely mimics the ac¬ 
tual spectroscopic LOWZ sample. The small differences can be ex¬ 
plained by small changes in the photometry used for target selection 
over the course of the BOSS survey. 


APPENDIX B: SDSS QUERIES 

B1 How to get photometric properties from DR8 

To obtain data from SDSS, we used the SDSS CasJobs tool 
First, we match objects from the redMaPPer catalog to SDSS data 
in RA and Dec using a 1.5 arcsecond match radius. The results of 
this match are stored in matchtable. We select the properties of 
the matched redMaPPer galaxies from SDSS DR8 using the SQL 
query: 

SELECT 

matchtable.search_id, matchtable.ra, 

matchtable.dec, p.modelMag_r, p.modelmag_i, 

p.modelmag_g, p.modelmag_z, p.modelmag_u, 

p.modeImagerr_u, p.modeImagerr_g, 

p.modelmagerr_r, p.modelmagerr_i, 

p.modelmagerr_z, p.cmodelMag_r, 

p.cmodelmag_i, p.cmodelmag_g, p.cmodelmag_z, 

p.cmodelmag_u, p.cmodelmagerr_u, 

p.cmodelmagerr_g, p.cmodelmagerr_r, 

p.cmodelmagerr_i, p.cmodelmagerr_z, 

p.extinction_r, p.extinction.!, 

p.extinction.g, p.extinction.z, 

p.extinction_u, p.petromag_r, p.psfmag_r, 

p.petror50_r, p.devrad.i, p.flags, p.objid 

INTO mydb.outputtable 

FROM mydb.matchtable 

JOIN PhotoPrimary as p on 

matchtable.matched_id=p.objid 

ORDER BY matchtable.search_id 


B2 How to get actual LRG target from SDSS 

To get the SDSS DR7 LRG target catalog (TargLRGs), we 
queried the SDSS DR7 database for all galaxies flagged as 
TARGE T_GALAXY_RE D: 

SELECT 

run, camCol, rerun, field, objID, ra, 
dec, modelmag.u, modelmag.g, modelmag.r, 

http://skyserver.sdss3.org/CasJobs/ 
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modelmag_i, modelmag_z, modelmagerr_u, 
modelmagerr_g, modelmagerr_r, modelmagerr_i, 
modelmagerr_z, extinction_u, extinction_g, 
extinction_r, extinction.!, extinction.z, 
petromag.r, psfmag.r, petrorSO.r, devrad.i 
into mydb. DR7_LRG_all 
from TargPhotoPrimary 
WHERE (primtarget & 

dbo . fPrimTarget ('TARGET_GALAXY_RED') ) != 0 
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